Observation of deviations from ideal gas thermodynamics in a trapped Bose-Einstein 

condensed gas 
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We have investigated experimentally the finite-temperature properties of a Bose-Einstein con- 
densed cloud of *^Rb atoms in a harmonic trap. Focusing primarily on condensed fraction and 
expansion energy, we measure unambiguous deviations from ideal-gas thermodynamics, and obtain 
good agreement with a Hartree-Fock description of the mixed cloud. Our results offer for the first 
time clear evidence of the mutual interaction between the condensed and thermal components. To 
probe the low-temperature region unaccessible to the usual time-of-flight technique, we use coherent 
Bragg scattering as a filtering technique for the condensate. This allows us to separate spatially 
the condensed and normal components in time of fiight, and to measure reliably temperatures as 
low as 0.2 r^" and thermal fractions as low as 10%. Finally, we observe evidence for the limitations 
of the usual image analysis procedure, pointing out to the need for a more elaborate model of the 
expansion of the mixed cloud. 

PACS numbers: 03.75.-b,03.75.Hh 
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I. INTRODUCTION 



Trapped, dilute atomic gases offer a new opportunity 
to study the interplay between quantum-statistical phe- 
nomena and interactions in Bose systems l^. A third 
ingredient, the external trapping potential, also plays a 
key role in understanding the properties of these gases 
At finite temperatures, it leads to spatially distinct 
condensed and thermal phases, a new behavior when 
compared to bulk quantum fluids, in which both compo- 
nents overlap everywhere. This spatial separation allows 
in particular a clear identification of the condensed frac- 
tion through absorption imaging, in stark contrast with 
superfluid *He, where condensed fraction measurements 
are only indirect 0. 

Much work has been devoted to the properties of con- 
densed gases at very low temperatures (much smaller 
than the critical temperature Tc), where the non- 
condensed fraction is negligible. Then, in the so-called 
Thomas- Fermi (TF) regime the static and dynamic 
behavior of the condensate is essentially determined by 
the interplay between the trapping potential and the 
atomic interactions. At higher temperatures, less than 
but comparable to Tc, a significant thermal component 
is also present, typically much more dilute than the con- 
densate. In this situation, the kinetic energy per thermal 
atom is larger than the mean-field energyjand deviations 
from ideal gas behavior are small 

iHillii Al- 

though several key thermodynamical properties, such as 
condensed fraction and average energy, are readily mea- 
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surable experimentally, a detailed comparison of experi- 
ments with finite-temperature theories of the interacting 
cloud is to our knowledge still lacking. 

The goal of this paper is to contribute to fill this gap, 
by investigating experimentally the thermodynamics of 
a trapped ^^Rb Bose gas below Tc, where both a con- 
densed and a thermal component are present. We fo- 
cus on two quantities, the condensed fraction and the 
expansion (kinetic plus mean-field) energy in the radial 
direction. We find that the condensed fraction Nq/N is 
significantly reduced with respect to the ideal-gas law, 
No/N = 1- {T/T^f, where Tc° is the ideal gas conden- 
sation temperature in the thermodynamic limit, and that 
the expansion energy is increased, even for the thermal 
atoms. 

Our measurements thus clearly exclude ideal-gas be- 
havior, and to assess the importance of interactions, we 
compare them with two mean-field theories of the in- 
teracting cloud that assume a condensate in the TF 
regime. The simplest one is the so-called "semi-ideal" 
model [lol ITU , which considers the thermal cloud as a 
quantum-saturated gas evolving in the combined trap- 
ping plus condensate-mean-field potential. Although 
the general trend of our observations is well reproduced 
by this model, finer details are not. To take interac- 
tions within the thermal cloud and between the thermal 
cloud and the condensate into account, we use a self- 
consistent Hartree-Fock (HF) description of the mixed 
cloud d El Q El El lii 111, which yields good 
agreement with the data. We are able to confirm exper- 
imentally its validity over a wide range of temperatures 
and atom numbers, and to show that despite the dilute- 
ness of the thermal cloud, its mean-field energy affects 
both the condensed and non-condensed atoms. Such sys- 
tematic measurements of the mutual interplay between 
the condensed and non-condensed components have not 
been reported before, although evidence for repulsion of 
the thermal atoms by the condensate has been provided 
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in [19j, through a careful study of the interface region 
btween the two components. 

Our measurements rely on the standard time-of flight 
technique to observe the mixed cloud. A limitation 
arises at low temperatures, where the condensate ap- 
pears broader than the thermal cloud in time of flight. 
This limits our measurements with this technique to 
T/T° > 0.3 (condensed fraction < 80 %). To overcome 
this limitation, we use coherent Bragg scattering to spa- 
tially separate the coherent and incoherent components 
of the cloud in time of flight. We apply it for low temper- 
ature thermometry, and measure in this way condensed 
fraction larger than 90 % and temperatures below 100 
nK (T/TO « 0.2). 

The paper is organized as follows. Section first de- 
scribes the experimental apparatus used for condensate 
production. Then we expose the widely used procedure 
of time-of-flight imaging, and introduce the improved 
method based on Bragg diffraction to filter the conden- 
sate out of thermal cloud. The key results of this paper 
are presented in sections IIIII and IIVI In section IIIII we 
present a measurement of the temperature dependance of 
the condensed fraction. In section Hvl the Bragg filter- 
ing scheme is applied to a measurement of the expansion 
energy of the thermal component. Both measurements 
clearly confirm the HF predictions in the whole temper- 
ature range. Further evidence for mutual interaction is 
provided in section where we analyze the behavior of 
the mixed cloud in time of flight and measure a compres- 
sion of the condensate axial length by the thermal cloud. 
For completeness, we review briefly in appendix ^ the 
three well-understood models to which we compare the 
data, namely the ideal gas model, the semi-ideal model 
and the self-consistent HF model. 



II. EXPERIMENTAL METHODS 

This section reviews the experimental techniques used 
in this work. After a brief description of our experimental 
apparatus Ijll A|) , we discuss in section III Bl condensed 
fraction and temperature measurements via time-of-flight 
absorption imaging, and point out to the limitations of 
this technique for high condensed fractions. In section 
III CI we show how they can be overcome using coherent 
Bragg scattering poll2l|. 

A. Condensate production 

Our experimental setup employs the standard combi- 
nation of laser and evaporative cooling to reach Bose- 
Einstein condensation in a sample of *^Rb atoms, spin- 
polarized in the |-F ~ l;mp — —1) hyperfine ground 
state. A Zeeman-slowed atomic beam loads a magneto- 
optical trap in 4 seconds. After spatial compression (125 
ms) and cooling in optical molasses (6 ms), the laser 
cooled sample is repumped in the \ F = 1) hyperfine 



manifold and loaded into an loffe-Pritchard type mag- 
netic trap. The trap, an iron-core electro-magnet, is 
a specific feature of our experiment [S^- The guiding 
of magnetic flux lines along the ferromagnetic body of 
the magnet produces strong (when compared to macro- 
scopic, coil-based devices) radial gradients (1.4 kG/cm), 
while using a modest current of 30 Amps. This pro- 
duces a cigar-shaped, very-elongated harmonic potential 
of the form Fext = MivHx"^ + y'^)/2 + MujIz'^/2. In 
this work, differently from |23l |. radial and axial trap- 
ping frequencies are respectively Lo^jlTx = 413(5) Hz and 
u:,j2-K = 8.69(2) Hz at a 10 G bias fleld. 

The combined compression and evaporation sequence 
typically lasts for 16 seconds, with the last part of 
the radio-frequency (rf) evaporation ramp considerably 
slowed down (to a ramp speed of 200kHz/s) and fol- 
lowed by a 1 s hold time in the presence of an rf shield. 
This is done in order to allow the system to cross slowly 
the quantum degeneracy threshold, and to relax towards 
equilibrium. In particular, non-equilibrium shape oscilla- 
tions that occur in such anisotropic traps upon conden- 
sation f23L 01 are strongly reduced. To ensure the re- 
producibility of the temperature in the experiment, the 
radio-frequency that empties the trap is measured ev- 
ery flve experimental cycles, and the rf shield frequency 
i^if is readjusted in real time to follow slow drifts of 
the bottom of the trap. In this way, the "trap depth" 
v-ci — vq is controlled within ±2 kHz. Since we measure 
^(y-sii — VQ)lk^T « 11 in this flnal evaporation stage, we 
estimate the temperature reproducibility to be ±10 nK. 



B. Time of flight imaging 

Informations on the atomic cloud are obtained through 
time-of-flight absorption imaging (see i25il for a detailed 
account of imaging techniques, and |2y] for details of 
the implementation in our experiment). This widely 
employed method will be referred to as the "standard 
method" in the following. Here, two important quan- 
tities are measured through a fit to absorption images, 
the condensed fraction and the temperature. To find the 
condensed fraction, we use the fitting method described 
in |25| , which assumes that the atomic density can be de- 
scribed by a parabolic TF profile describing the conden- 
sate plus an ideal Bose-Einstein distribution with zero 
chemical potential describing the thermal cloud. The 
condensed number A'o is taken to be the number of atoms 
integrated under the parabolic profile, while the total 
atom number N is found by integration over the entire 
profile. 

The temperature is measured through an independent 
fit to a thermal Bose-Einstein distribution, restricted to 
the wings of the thermal distribution only, in order to 
minimize the effect of interactions (see '3, and section 
II V|) . One typically assumes that the high-energy atoms 
selected by this procedure behave as if the gas were ideal, 
and extracts effective temperatures and from the 
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cloud sizes, for instance kuTx = Mujj_R^y^/[l + {uj±t)'^] 
with i?th the radial size. We observe a systematic vari- 
ation of the temperature Tj measured along the long 
axis with the size of the exclusion region. This effect 
was more pronounced with increasing condensed fraction, 
and disappeared above Tc- On the contrary, the radial 
temperature was barely affected by the actual size of 
the exclusion region, provided it was chosen larger than 
the condensate radius and sufficiently small to conserve 
a reasonable signal-to- noise ratio (5% or less variation). 
For this reason, we infer the initial temperature from the 
radial value only |28| . This thermometry procedure 
assumes a nearly-ballistic expansion. CoUisional effects 
in time-of-flight may invalidate this assumption, and it is 
important to estimate their importance to quantify the 
accuracy of our measurements. We defer this discussion 
to sections HVl and IVl 



C. Selective displacement of the condensate using 
coherent Bragg scattering 

The method described in the last subsection rests on 
the clean distinction between the condensed and thermal 
components. However, for cold samples (fceT < fj,), the 
condensate radius after time of flight is larger than the 
extension of the thermal cloud. The details of the ther- 
mal distribution are therefore "buried" under the conden- 
sate density profile, and a reliable fit is not possible. We 
have found no difficuhies for No/N - 60% (T/T^ - 0.5), 
while still being able to extract valuable informations 
for condensed fractions close to 80 % (T > 0.3 Too). 
For higher condensed fractions, the signal-to- noise on the 
wings is too low to perform a reliable fit. 

In this section, we describe an improvement of the 
time-of-flight technique that exploits the dramatic differ- 
ence between the two components in momentum space 
0, l20l l2lj | . Our primary motivation for introducing this 
method is to improve the resolution at low temperatures, 
a key advantage in the condensed fraction measurements 
presented below. The principle of the technique is to 
transfer the whole condensate to a center-of-mass mo- 
mentum state with two photon recoil velocities (approx- 
imately 1.1 cm/s in our case), while leaving the thermal 
atoms essentially unaffected. As shown in Fig. ^ after 
a subsequent time of flight (typically 20 ms), the con- 
densate has moved away from the center of the thermal 
cloud, allowing us to perform separate fits and revealing 
previously hidden information about the thermal com- 
ponent. This technique is related to the filtering scheme 
used in [S^ to study collective modes in the presence of 
a weak optical lattice at finite temperatures. 

Momentum transfer is achieved by exposing the sam- 
ple to two counter-propagating laser beams detuned with 
respect to each other. The moving lattice formed this 
way can Bragg diffract the atoms, promoting them to a 
two-recoil momentum state. In principle, Bragg scat- 
tering is sensitive to the atomic momentum distribu- 



tion [23, [^. However, if the spectral width of the 
Bragg pulse, Fourier-limited by the pulse duration T to 
Ai'p '-^ 1/T, is much larger than the condensate Doppler 
width Avq ^ 2ur/L, the momentum sensitivity of Bragg 
diffraction is lost and almost complete transfer of the 
condensate (tt pulse) is achieved. On the other hand, if 
the thermal component spectral width, Avth ^ 2wr/At 
(At = \J 2Tih? / mk-oT is the thermal De Broglie wave- 
length), is much larger than At'p, Bragg diffraction is 
then a momentum-selective process, and most of the ther- 
mal atoms, being off-resonant with the Bragg beams, 
stay at rest in the laboratory frame. This filtering tech- 
nique thus makes possible independent manipulation of 
the condensed and non-condensed atoms. 

The experimental setup that we use to generate the 
Bragg beams is identical to the one used in [SJ. The 
TT Bragg pulse is applied after sudden switch-off of the 
trapping potential and a 2 ms time of flight, to avoid 
s-wave collisions between atoms with different momenta 
|30j | . The Bragg beams are counter-propagating and par- 
allel to the elongated axis of the trap. Since the conden- 
sate mean-field energy is then almost entirely released in 
the transverse directions, perpendicular to the diffraction 
axis, the ratio between the condensate and thermal cloud 
spectral widths is 

and remains smaller than one even at very low temper- 
atures. Note that this would not have been the case for 
Bragg beams parallel to a radial axis, where the momen- 
tum width is larger by a factor [uj^/uzf due to the the 
released mean field energy. 

To find the condensate fraction, we simply count the 
number of diffracted atoms and identify it with the con- 
densed number. To obtain a meaningful measurement of 
the condensed fraction, two conditions have to be met. 
First, almost all condensate atoms should be diffracted. 
Transfer efficiencies as high as 97 % were observed for the 
coldest clouds we have produced, indicating that only a 
few percent of the condensate population remains at rest; 
this is comparable to the ±2 % root-mean-square fluctu- 
ations of the diffraction efficiency that we have measured 
experimentally. Second, the diffracted number of thermal 
atoms should be small to obtain a sensible measurement 
of the condensed number. We estimate the diffracted 
fraction of thermal atoms /th in the impulse approxima- 
tion discussed in |3lj |. In this approximation, well justi- 
fied for the low densities considered here, the number of 
atoms diffracted from the thermal cloud is taken to be 
proportional to the thermal part of the momentum dis- 
tribution, integrated over the directions perpendicular to 
the axis of the Bragg beams. Treating the thermal cloud 
as an ideal Bose gas in a first approximation, we find for 
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a resonant tt pulse 

, _ TT^ [35/2(1) - g5/2 {exp [-16(i/R/Azyth)^] }] 
■^"^ ~ AAm^T 53/2(1) 

(2) 

Here and below, the Bose functions are defined by 
ga{u) — X]j>i '""'/i"- The second term in brackets 
arises from off-resonant excitations, which are negligi- 
ble for the condensate but sizeable for the thermal com- 
ponent. The quantity fth varies from 6 % at 400 nK 
(T - 0.8 TO, or Nq/N - 20 %) to 20 % at 150 nK 
(T - 0.2 r°, or No/N - 90 %). The number of thermal 
atoms in the diffracted peak is given by (1 — iVo/A^)/th- 
All in all, we estimate from this simplified calculation 
that the Bragg filter is a useful technique in a window 
5% < Nq/N < 95%, the upper bound being set by the 
first condition (complete condensate transfer) and the 
lower by the second (number of condensed atoms larger 
than number of thermal atoms in the diffracted peak). 
Outside of this window, it becomes suspect to identify 
unambiguously the diffracted order with the condensate. 
Despite of this, the technique represents a significant im- 
provement over the standard method, as we will sec in 
the next sections. 



III. CONDENSED FRACTION OF THE 
INTERACTING GAS VERSUS TEMPERATURE 

With the tools of section |n] in hand, we can investi- 
gate the behavior of the condensed fraction as a function 
of temperature. This quantity is of primary importance: 
the possibility to measure it directly in trapped gases is 
in stark contrast with the situation in traditional, bulk 
superfluids where such measurements are intrisically dif- 
ficult Q. We have done this measurement in two steps, 
first by using the standard method and second by tak- 
ing advantage of the enhanced resolution of the Bragg 
filtering scheme. 

A. Standard time-of flight measurements 



increase the condensed fraction to diminish the exchange 
interaction energy among excited states. 

An interesting property, first pointed out in j^^, is 
the scaling behavior of all thermodynamical quantities, 
which depend only on the reduced temperature T/T^ and 
on the parameter 77 which controls the magnitude of two- 
body repulsion .411]. 

Alternatively, one can express 77 as the ratio between the 
two characteristic lengths in the uniform problem, the 
scattering length a and the De Broglie wavelength Aq for 
T = T^, as 7? « 1.07(a/Ao)2/^ The power 2/5 reflects 
the presence of the trapping potential. A typical value 
in our experiment is ry = 0.5, while the experiment of 
(for instance) corresponds to ry 0.3 — 0.4. Thus, the 
effect of two-body interactions is stronger in the work 
reported here, which explains to some extent the clarity 
with which we observe deviations from ideal gas behavior. 

In the experiment, because the total number N drops 
with T due to evaporation and losses, rj slightly decreases 
across the temperature range, from 0.51 above to 0.47 
at low temperature (see Fig. 12; and d). To compare with 
theory, we use the average number of atoms across the 
data set, N « 1.2 x 10^, and the corresponding rj = 0.49. 
This docs not lead to a discernible variation of the HF 
prediction at the scale of the graph, because of the very 
weak N^^^^ dependence of rj. This behavior emphasizes 
the scaling behavior exhibited by trapped Bose gases with 
large atom numbers 33]. 

Finally, we note that the data in the vicinity of T^, 
shown in Fig. emphasize the necessity of the full HF 
treatment of the thermal component to understand quan- 
titatively the thermodynamic properties. As a matter of 
fact, the semi-ideal model (dashed line) predicts a con- 
densed fraction systematically higher than the one we 
observe. This is clear evidence for mutual interaction 
between the condensed and non-condensed components. 
We will return later to this point, which is an important 
conclusion to be drawn from this work. 



Using the standard analysis technique (i.e. no Bragg 
filtering), we investigate in this section a temperature 
interval ranging from T « 1.1 T° down to T « 0.3 T°, 
with approximately 2 x 10^ atoms at the transition. Fig.|5] 
shows the measured condensed fraction as a function of 
T/T^. Each point results from an averaging over several 
(typically five) realizations under identical conditions. 
The expectation for an ideal gas, including finite-size 
effects lies distinctly above our experimental data. 
The difference can be attributed to interactions, as shown 
by the far better agreement with the self-consistent HF 
calculation (solid line). The observed reduction of the 
condensed fraction contrasts with the homogeneous case 
|32j ]. where the condensed and normal component over- 
lap everywhere, and where it is energetically favorable to 



B. Enhanced resolution of low thermal fraction 

As discussed earlier, at low temperatures, T < 0.4Tc, 
the standard procedure is unable to extract faithfully the 
properties of the thermal cloud. Thanks to the Bragg 
filtering technique introduced in section III CI this diffi- 
culty can be overcome and very small thermal fraction 
can be detected. The condensed fraction measured this 
way is plotted in Fig. ^jp, along with the rj = 0.47 curve 
that corresponds the value N ^ 8 x 10^ for this set of 
measurements. Again, we find good agreement with the 
prediction of the HF model within our uncertainty, even 
at very low temperatures. One could wonder whether the 
contribution from collective excitations (quantum and 
low-energy thermal depletion) could be measured by this 
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technique (the analogous of the phonon regime in super- 
fluid ^He). Unfortunately, according to the estimations 
of , they are always small when compared to the con- 
tribution of single-particle excitations in the temperature 
range wc explore, and compared to the estimated sensi- 
tivity of the Bragg technique fsection III Only by in- 
creasing the diluteness parameter \J no(0)a'^ significantly 
could this regime become observable experimentally with 
the techniques described here. 



IV. MEASUREMENT OF THE OVERALL AND 
THERMAL EXPANSION ENERGIES 

Another quantity that can be measured from time-of- 
flight expansion is the release energy Q, H, |^ , the sum 
of the kinetic and interaction energy released at the trap 
cut-off and available for the expansion of the whole cloud. 
In an anisotropic trap such as ours (lo^ju)^ ~ 51), almost 
all the interaction energy converts into radial expansion 
velocity. The radial expansion of the cloud for t 3> ■, 
observed in the y direction, proceeds at an overall speed 
fy, fixed by the expansion energy 

i?, = iM«2 = ii?kin + ^i?i„t, (4) 

which, if scaled by the characteristic Nk^T^, is a univer- 
sal function of rj and T/T^. We discuss in Appendix IXI 
how to calculate the kinetic and interaction energies in 
the HF approximation. 

Experimentally, one measures Ey directly from the 
root mean square cloud radius, according to {y^) = Vyt 
0, 15 , without resorting to a detailed fitting model. This 
expression assumes negligible relaxation between the ax- 
ial and radial degrees of freedom. We have plotted the 
measured value of Ey in Fig. for same data as those 
shown in Fig. \3p, where the Bragg filter has been used 
[r] K, 0.47 for these data). As expected, the kinetic energy 
of the thermal cloud dominates close to Tc, with a small 
contribution of the mean-field energy, whereas the inter- 
action energy of the condensate E^ is the most important 
term at low temperature. The HF curve connects these 
two limiting cases, and reproduces well our observations. 

A key advantage in our situation is the ability to an- 
alyze separately the condensate and the thermal cloud 
thanks to the Bragg filtering scheme, and therefore to 
measure the release energy of the thermal cloud alone. 
To avoid the condensate, we estimate the thermal cloud 
rms radius from a fit to a radial cut to the image (see in- 
set in Fig. Od). The release energy of the thermal cloud 
measured this way is shown in Fig. (SJ^, together with the 
calculated value, 

i?r^ =^i^ki„+^(Co+cth)- (5) 

The quantities Cth represent the scaled mean- field en- 
ergy corresponding to the repulsion felt by a thermally 



excited atom due the condensate and to the remaining 
thermal atoms, respectively. They are defined more pre- 
cisely in Appendix El Although the difference with the 
non-interacting curve is less pronounced at very low tem- 
peratures, where the thermal energy is very low, close to 
Tc, these results emphasize the important role of inter- 
actions and the good agreement with HF theory once 
again. 



V. FURTHER EVIDENCE FOR MUTUAL 
INTERACTION BETWEEN THE CONDENSED 
AND THERMAL COMPONENTS 

From images taken employing the Bragg filtering 
scheme, it is also possible to examine a radial cut to 
the profile of the thermal cloud, as shown in Fig. Ql^- 
As the initial kinetic energy of the thermal cloud is typi- 
cally much larger than its mean-field energy, one expects 
that the density distribution after time of fiight reflects at 
least approximately the initial momentum distribution, 
which does not display the "hole" present in the density 
distribution in the trap (see Appendix 1X1 for further dis- 
cussion). This is indeed the case: the measured profiles 
show a monotonic behavior near the center of the cloud 
p[^. Nevertheless, the density distribution we observe is 
somewhat flatter in the central region than the ideal gas 
distribution used in the analysis, as shown by the resid- 
uals of a flt to the radial proflle (Fig. ^Jd). Above Tc, 
the effect disappears, which indicates that a flatter pro- 
file is not simply an artifact of our measurement method. 
Furthermore, such a behavior is to be expected if the con- 
densate mean field repels the thermal cloud in the early 
stage of the expansion, since interactions tend in general 
to make the density profile more uniform. 



A. Evidence for non-ballistic expansion 

As stated in section^ this repulsion effect is not taken 
into account in our fitting procedure, which assumes an 
ideal Bose-distribution to fit the profile and ballistic ex- 
pansion to deduce temperature from the cloud sizes. To 
investigate further the validity of the analysis, we begin 
by plotting in Fig. and the aspect ratio of both 
components, as a function of the reduced temperature. 
One sees from these graphs that the simple model of a 
TF condensate on top of an ideal thermal background 
is not sufficient to account for the data. Indeed, for our 
trapping frequencies and for the time of fiight t = 22.3 
ms used in these measurements, one would expect from 
this model an aspect ratio of 1.17 for the condensate and 
0.77 for the thermal cloud. Both deviate from these val- 
ues and vary with temperature, indicating that the ex- 
pansion dynamics is more complex than assumed by the 
analysis model. Note that although this model fails to 
describe fully the expansion dynamics, the observed de- 
viation from ballistic expansion remains small. 
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As already pointed out, the obvious weakness of the 
analysis model is the neglect of collisional effects. We 
recall that the expansion model is basically motivated 
by the absence of a more elaborate theory to which wc 
could compare our observations. It is however of interest 
to quantify the error level on temperature measurements, 
which we will do in the remainder of this section at the 
gross estimate level. Interaction-driven forces that affect 
the expansion can be divided into two distinct classes 
|34| : hydrodynamic forces on one hand, predominant 
above Tc in |27j . and mean- field repulsion on the other, 
which play a minor role above Tc, but becomes increas- 
ingly important with decreasing T due to the presence 
of the condensate Both effects are roughly speak- 

ing comparable in magnitude, and from our observations 
above and HF calculations (see section^, we estimate 
an upper bound on the systematic error on temperature 
measurement, from 10 % close to Tc (mostly due to hy- 
drodynamic behavior) to 20 % well below (mostly due 
to repulsion by the condensate), compatible with the ob- 
served deviation from the ideal gas expansion. 

This is an upper bound because the temperature is 
found through a fit to the wings of the distribution the 
only, which is expected to reduce effect of mean-field re- 
pulsion. The close agreement we find between the cal- 
culated release energy and the measured one in Fig. |3|d, 
where we recall that the radial size used to deduce the ex- 
pansion energy was found through a fit to the full profile, 
seems to indicate that fitting to the wings avoids count- 
ing most of the repulsion energy in the determination of 
the temperature, as expected. However, this need not 
be true for the hydrodynamic relaxation, which affects 
directly the momentum distribution. Just as above Tc, 
hydrodynamic forces are thus likely to be the dominant 
source of systematic error in most of the temperature 
range in elongated traps such as ours. 



B. Compression of the condensate by the thermal 
cloud 

We conclude this section by examining the axial length 
of the condensate, which is reasonably immune to the ef- 
fects discussed above because the axial expansion of the 
condensate is very slow: the measured length thus stays 
close to the in-trap length. From condensed fraction mea- 
surements, it is clear that the semi-ideal model is not suf- 
ficient to reproduce our results, meaning that the mutual 
interaction between condensed and non-condensed atoms 
is observable. This also seen from the length of the con- 
densate inferred from the two-component fit, which is re- 
duced when compared to the TF length calculated with 
the number of condensate atoms we measure. A quantita- 
tive comparison can be made only by taking into account 
the slow axial expansion. In the absence of a complete 
theory, we assume that the axial length is rescaled from 
the equilibrium length by the same factor as a conden- 
sate in the TF regime, &2 « 1 +'KU}1t/2u)± « 1.04 for our 



parameters |35l ISq . With this assumption, we find our 
data to be in reasonable agreement with the HF length, 
whereas the TF prediction is found systematically too 
high (see Fig. EJ. This reduction can be attributed to 
the thermal cloud compression discussed in appendix 1X1 
at equilibrium, the shell of thermal atoms surrounding 
the condensate exerts a force towards the trap center, 
reducing its extension when compared to a "free" , TF 
condensate. 



VI. CONCLUSION 

In this paper, we have investigated experimentally the 
thermodynamics of a trapped, interacting Bose gas over 
a wide range of temperatures, from down to 0.2 T^. 
We have used the standard time-of-flight analysis, com- 
plemented by the use of coherent Bragg scattering to fil- 
ter the condensate out of the thermal cloud. The lat- 
ter technique allows us to reach lower temperatures and 
higher condensed fractions than those accessible by the 
usual method. We have investigated primarily two quan- 
tities, the condensed fraction and the (radial) release en- 
ergy. The data display without ambiguity an interact- 
ing gas behavior, and are in agreement at a few percent 
level with a Hartree-Fock description of the mixed cloud. 
This highlights the pertinence of this simple description 
of the interacting, trapped Bose gas at finite tempera- 
ture. Moreover, this gives evidence for mutual interac- 
tion between the condensate and the thermal cloud close 
to Tc. Although these effects are small, they are measur- 
able and should be taken into account in precise compar- 
isons to the theory of finite temperature, Bose-Einstein 
condensed gases. 

A more thorough quantitative test of theory for ground 
state occupation and similar thermodynamic properties 
is however hindered by the lack of theory to understand 
the expansion of a mixed cloud. For instance, the data 
in Fig. lie systematically slightly below the theory, 
which may be due to an underestimation of the con- 
densed fraction or the temperature: as described above, 
we use an ideal Bose distribution to fit the thermal com- 
ponent. This assumption contradicts the principal con- 
clusion of this paper, that ideal gas approximations are 
far from sufficient at the desired level of accuracy. Al- 
though the experiments described in this paper suggest 
an accuracy on temperature of order 10 % at least, we 
have also shown that the density profile deviates from 
an ideal gas near the center of the cloud, and stressed 
that hydrodynamic and mean-field effects in the expan- 
sion are not properly accounted for. Systematic errors 
may thus still be present, and a better accuracy is not 
guaranteed with the methods used in this paper. Nu- 
merical work, for instance along the lines of |33|, niay 
help study the expansion of mixed clouds and improve 
the standard analysis procedure. 
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APPENDIX A: SIMPLE MODELS OF A 
TRAPPED BOSE-EINSTEIN CONDENSED GAS 
AT FINITE TEMPERATURES 

In this Appendix, we briefly summarize the three dif- 
ferent models to which we compare our observations. 
We use as a temperature scale the critical temperature 
of an ideal gas in the thermodynamic limit's!, ^bT^ = 
huj[N / C,{i)Y^'^ , where C is the Riemann function, and 

—2/31/3. 

Lu = uj ^ uJz is the geometrical mean of the trapping fre- 
quencies. We make three key assumptions to simplify the 
theoretical description. For the condensate, we suppose 
that the condition Nqg/W 3> 1 holds, where W — y^h/Mu) 
is the mean ground state width and a the scattering 
length. This ensures that the condensate is in the TF 
regime. At very low temperatures, where the condensed 
fraction is almost unity, this gives the condensate density 
no = riTF as 

"TF(r) = [/iTF-Kxt(r)]/C/, (Al) 

where U = Anh^a/M is the mean- field coupling con- 
stant. The TF chemical potential is /^tf = (fi^/2) x 
{15Noa/a)^/^. For the thermal cloud, we assume ^ 
huj± for the semi-classical approximation to hold, and use 
a mean-field description which treats thermal atoms as 
independent particles evolving in a self-consistent static 
potential Kff(r) (see below). The distribution function 
in phase space for thermal particles then reads 

fi^^P)^ ^0iniJ)-,]-r ^^^^ 

with the semi-classical hamiltonian 7i(r, p) = p^/2Af -|- 
Vcs{r), the chemical potential ^, and /3 = l/ksT. The 
particle density distribution for the thermal component, 
nth, is found by integegration over momenta and reads 

«th(r) = 53/2 {exp [(3 (^ - V^ff (r))]} . (A3) 

Each model detailed below is thus specified by the precise 
form of Vcs'- in the ideal gas model, Vcs reduces to the 
trapping potential, in the semi-ideal model, it includes 
the mean-field of the condensate only, and in the HF 
model, it also takes the mean field of the thermal atoms 
into account. 



1. Ideal gas model 

The simplest approximation neglects all interactive 
contributions to the effective potential, which reduces to 
the trapping field, 

Kff - M = Kxt - (A4) 

The thermodynamic quantities follow the ideal-gas laws 
(see for instance j^), and impose in particular /i = 

below Tc and ^ = 1 - (^^) for T < T^. This ideal 
gas description is thus incompatible with the existence 
of a TF condensate, that imply fi = /itf > 0, and it 
should be considered as a reasonable approximation only 
for k^T ^ jjL. 



1. Semi-ideal model 

The repulsion of the thermal cloud by condensed atoms 
is taken into account in the so-called "semi-ideal" model 
[lTI| . which considers a TF condensate containing 
A^o < ^ atoms plus a quantum-saturated, ideal thermal 
gas moving in the combined trapping plus condensate 
mean field potential, 

yoff(r) - M = Vcxt(r) + 2[/nTF(r) - Mtf- 

The factor of 2 accounts for exchange collisions between 
atoms in different quantum states [3| . Note that the con- 
densed atom number should be found self-consistently 
for a given T, under the constraint that the total atom 
number is fixed. 

The semi-ideal model correctly predicts the most im- 
portant feature in the static density profile of trapped in- 
teracting Bose gases. As soon as the condensed fraction 
is larger than a few percent, the density of the conden- 
sate greatly exceeds the density of the thermal compo- 
nent. Therefore, the condensate mean field is stronger 
and repels the thermal cloud from the center of the trap, 
digging a hole in the thermal density distribution. 



3. Self-consistent Hartree-Fock model 

Although qualitatively correct and appealing because 
of its simplicity, the semi-ideal model is not sufficient to 
describe precisely our experiments. Interactions among 
thermal atoms and the backaction of the thermal cloud 
on the condensate have to be taken into account, which 
we do here in the HF approximation [l3 . Il3 . \\A . llSl lig . 
IT^ ll^ , corresponding to a self-consistent potential 

Kff(r) - M = Kxt(r) + 2C/7io(r) + 2C/7ith(r) - y.. (A5) 

The equilibrium condensate density is no longer simply 
given by the TF profile, but depends also on the thermal 
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density through 



no[r) 



U 



(A6) 



A self-consistent numerical solution with a fixed atom 
number N fixes the finite T chemical potential, fi = 
Uno{0)+2UnthiO), from which all other quantities can be 
determined. For instance, the interaction energy, which 
contains mixed terms describing the mutual influence of 
the condensate and the thermal cloud, can be written as 
Eint =Eo + fcBr(2Co + Cth), where 



Eq 

Co 

Cth 



u 



u 

U 



d(3)r noirf, (A7) 
d(3)r no(r)nth(r), (A8) 
d(3)r (nth(r))' . (A9) 



The kinetic energy is entirely due to the thermal cloud 
in the TF approximation, and reads 



l-d'^^^rS^^p — f(r,p). (AfO) 



Numerical solution of the HF model displays two addi- 
tional features compared to the semi-ideal model. First, 
mean-field interactions lower the critical temperature for 
Bose-Einstein condensation p?! Is^ . Second, there is a 
back-action of the thermal cloud on the condensate: the 
mean-field exerted by the shell of thermal atoms sur- 
rounding the condensate acts in return to compress it, 
increasing its density and reducing its axial length Lq 
according to 



25 



'^o(O) + ^ ( 53/2 



g"0(O) 



53/2 [1]) }( 



All) 



This compression effect is directly observed in section 
IV Al and indirectly through the measurements of con- 
densed fraction in section imi 

A further step forward would be to include collective 
effects in the model. However, reference ,38] points out 
that low-energy, collective excitations cause a minute 
change in the thermodynamic properties of the system 
even at relatively low temperatures, T < fi. Another 
app roach, based on Quantum Monte-Carlo calculations 
[33. l40l | , has confirmed that the HF approximation could 
reproduce the thermodynamics of the trapped clouds to 
a very good accuracy. 
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FIG. 1: Bragg difFraction as a condensate filter. Figures (a), 
(b) and (c)correspond to 92, 85 and 30 % condensed fraction, 
respectively {T/T° ^0.2, 0.35 and 0.8). Two-dimensional ab- 
sorption images are shown, with a cut along the direction of 
the trap weakest axis ( "axial cut" ) . The top images show reg- 
ular absorption images for a 24.27 ms time of flight (dashed 
lines in the cut) . The bottom images, also shown as the solid 
line in the axial cut, corresponding to the same temperatures 
and atom numbers within experimental reproducibility, have 
been taken after applying a moving optical lattice tuned to 
realize a vr Bragg pulse, transferring two photon recoils to al- 
most all condensed atoms (velocity 1.1 cm/s). The distance 
travelled by the condensate is verified to be 250 ^m, corre- 
sponding to a free flight of 22.27 ms after the Bragg pulse. As 
can be seen, the thermal cloud is barely affected, due to its 
much larger extent in momentum space. 
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FIG. 2: Condensed fraction as a function of reduced temper- 
ature. (a),(b) Experimental data after averaging are shown 
as filled circles, with statistical error bars. Data in (a) were 
taken using the standard time-of-flight technique, while data 
in (b) were measured with the Bragg filtering scheme and ex- 
tend to lower temperature, down to T « 0.2 Tco whereas the 
standard method is limited to T « 0.4 Tco. Lines show theo- 
retical expectations according to an ideal gas calculation, in- 
cluding finite-size effects (long dashed), a "semi-ideal" model 
that neglects interactions within the thermal cloud (short 
dashed), and a self-consistent HF calculation (solid). Due to 
a different average atom number, the parameter rj that con- 
trols the importance of interactions is different in each case, 
fj = 0.49 and rj — 0.47 for (a) and for (b), respectively. The 
variation in total number due to evaporative cooling across 
the data set shown in (a) is shown in (c), and the corre- 
sponding variations of rj in (d) (Note the vertical scale, ex- 
tending over no more than 5 % of the average value). The 
average rj = 0.49 is shown by the dashed line. In (e), we 
show an enlargement of (a) around the critical temperature, 
to highlight the importance of making the full HF calculation 
to reproduce the trend seen in the data. 




FIG. 3: Expansion energy in the radial x direction as a func- 
tion of temperature. Filled circles in (a) correspond to the 
whole cloud, and open circles in (b) to the thermal compo- 
nent only. The data are taken from the same set as in Fig. 
The solid lines on the graphs show the same quantities pre- 
dicted by the self-consistent HF model, with rj = 0.47. The 
dotted line is the expansion energy (kinetic only) of an ideal 
cloud. 
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FIG. 4: Radial density profile of the thermal cloud. The 
dashed line in the inset shows the axis along which the profile 
are taken, (a) Radial cut of the column density profile of 
the thermal cloud after 24.3 ms time of flight (solid) and the 
best fit to an ideal Bose-Einstein ditribution (dashed). The 
Bragg filter has been employed to separate the thermal and 
condensed components. The condensed fraction is indicated 
in each case, (b) Residual of the fit for each case in (a). 




FIG. 5: Deviation from ballistic expansion. The aspect ratio 
of the condensed (a) and thermal (b) components of mixed 
clouds, after 22.3 ms of free expansion, are plotted as a func- 
tion of the reduced temperature. The aspect ratio for non- 
condensed clouds, analyzed in more details in |2^. are also 
shown for comparison. The horizontal dashed lines indicate 
the aspect ratio of a TF condensate and an ideal thermal gas, 
respectively, and the vertical dotted lines show the critical 
temperature including mean-field and finite-size effects. 
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FIG. 6: Compression of the condensate by the thermal cloud. 
The axial length of the condensate as measured in the absorp- 
tion images is shown, with the TF (dashed) and the HF pre- 
dictions for the trapped condensate (dotted) and after rescal- 
ing by the same factor as a TF condensate (solid). 



